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ABSTRACT 

We present Spitzer/IRAC photometry on the star-forming H II region N 90, related to the young stellar 
association NGC 602 in the Small Magellanic Cloud. Our photometry revealed bright mid-infrared sources, 
which we classify with the use of a scheme based on templates and models of red sources in the Milky Way, 
and criteria recently developed from the Spitzer Survey of the SMC (Bolatto et al. 2007) for the selection 
of candidate Young Stellar Objects (YSOs). We detected 57 sources in all four IRAC channels in a 6.'2 x 
4f 8 field-of-view centered on N 90; 22 of these sources are classified as candidate YSOs. We compare the 
locations of these objects with the position of optical sources recently found in the same region with high- 
resolution HST/ACS imaging of NGC 602 by Schmalzl et al. (2007), and we find that 17 candidate YSOs 
have one or more optical counterparts. All of these optical sources are identified as pre-main sequence stars, 
indicating, thus, ongoing clustered star formation events in the region. The positions of the detected YSOs 
and their related PMS clusters give a clear picture of the current star formation in N 90, according to which 
the young stellar association photo-ionizes the surrounding interstellar medium, revealing the H II nebula, and 
triggering sequential star formation events mainly along the eastern and southern rims of the formed cavity of 
the parental molecular cloud. 

Subject headings: Magellanic Clouds — HII regions — open clusters and associations: individual (NGC 602) 
— stars: formation — stars: pre-main sequence — Hertzsprung-Russell diagram 



1. INTRODUCTION 

The Spitzer Space Telescope opens a unique mid- and 
far-infrared window to the exploration of massive star for- 
mation not only in our galaxy but also in our neighbor- 
ing dwarf galaxies, the Large and Small Magellanic Cloud 
(LMC, SMC). Stellar associations in both the Magellanic 
Clouds (MCs) contain the richest sample of young bright stars 
in these galaxies, and consequently our knowledge on their 
young massive stars has been collected from studies of such 
stellar systems (see e.g. Massey 2002). Furthermore, al- 
most every MCs association coincide with one or more H II 
regions, as they have been earlier cataloged (Henize 1956; 
Davies et al. 1976), and it has been shown that the measured 
Ha fluxes are in excellent agreement with those expected 
from the ionizing flux of the detected massive stars (Massey 
1993). However, MCs associations are not mere aggregates of 
young bright stars alone, but they also host large numbers of 
faint pre-main sequence (PMS) stars, as recent HST studies 
showed for both the LMC (Gouliermis et al. 2006) and the 
SMC (Nota et al. 2006). Both, the existence of H II regions 
and PMS stars in stellar associations of the MCs indicate that 
star formation may be still active in their vicinity, and Spitzer 
offers a unique opportunity to test this hypothesis. 

Indeed, using Spitzer observations, Jones et al. (2005) re- 
ported four Young Stellar Objects (YSOs) in the H II com- 
plex N 159 (located in the vicinity of the association LH 105; 
Lucke & Hodge 1970) in the LMC. Chu et al. (2005) obtained 
Spitzer observations of the super-bubble N 5 ID (associations 
LH 51/54) in the LMC and they found three YSOs projected 
within the super-bubble. Images from the Hubble Space Tele- 
scope (HST) in Ha and [S II] allowed them to identify dust 
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globules associated with two of them. The third YSO was 
associated with the first Herbig-Haro object detected outside 
the Galaxy. Furthermore, Meixner et al. (2006) in their pre- 
sentation of the first results from the Spitzer Space Telescope 
Survey of the LMC Legacy Project Surveying the Agents of a 
Galaxy's Evolution (SAGE), show as a test-case epoch 1 data 
of the area of two H II regions, N 79/N 83 (associations LH 
2/LH 5). They found that MIPS 70 ^m and 160 pm images of 
the diffuse dust emission of the regions reveal a distribution 
similar to the gas emissions, especially the H I 21 cm emis- 
sion, and they note that the measured point-source sensitivity 
for the epoch 1 data is consistent with expectations for the 
survey. 

Bolatto et al. (2007), recently, presented their results from 
the Spitzer Survey of the Small Magellanic Cloud (S 3 MC), 
which imaged the star-forming body of the SMC in all seven 
MIPS and IRAC wave bands. These authors compiled a 
photometric catalog of 400,000 mid- and far-infrared point 
sources in the SMC, from which they identified 282 bright 
candidate YSOs as bright 5.8 /im sources with very faint opti- 
cal counterparts and very red mid-infrared colors ([5.8]-[8.0] 
> 1 .2 mag). More recently, Simon et al. (2007) used the same 
observations and applied spectral energy distribution (SED) 
fits based on models by Whitney et al. (2003a, b, 2004) to 
enhance the catalog of candidate YSOs in the brightest H II 
region in the SMC, N 66. 

The present paper deals with mid-infrared imaging from 
the Spitzer Space Telescope of the H II region LHA 1 15-N90 
(Henize 1956) or in short N 90, located in the wing of the 
SMC. This region, also cataloged by Davies et al. (1976) 
and named DEM S 166, is related to the young star cluster 
NGC 602, which was classified as a stellar association by 
Hodge (1985; system No 68 in his catalog). We present the re- 
sults of our investigation on the candidate YSOs in the region 
of NGC 602/N 90, identified with our photometry according 
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FIG. 1 . — Spitzer image of the region surrounding NGC 602/N 90. IRAC three-color image is shown, with IRAC-1 3.6 fj,m in blue, IRAC-2 4.5 fim in green, 
and IRAC-4 8 fim in red. 



to criteria based on galactic templates and recent Spitzer ob- 
servations of other star forming regions in the SMC. 

Specifically, in §2 we discuss the retrieval of the Spitzer 
data used in our study and their photometry. General descrip- 
tion of the mid-infrared characteristics of the region is given 
in §3, and in §4 the point sources detected with our pho- 
tometry, and the selection of candidate YSOs are presented. 
The recent release of high-resolution data on NGC 602/N 90 
obtained with the Advanced Camera for Surveys on Hubble 
Space Telescope allowed us a detailed photometric study of 
the stellar association and its surroundings down to the faint 
red pre-main sequence stars of the system (Schmalzl et al. 
2007). In §4 we also combine the results from both Hubble 
and Spitzer photometries to identify the optical counterparts 
of the candidate YSOs found in the region. Based on our find- 
ings we discuss the current star forming process that probably 
took place in the region of NGC 602/N 90 with emphasis to 
the clustered behavior of star formation in §5. Finally, general 
conclusions of this study are given in §6. 

2. DATA RETRIEVAL AND PHOTOMETRY 

The imaging data of our study have been obtained with the 
Infrared Array Camera (IRAC; Fazio et al. 2004) on-board 
Spitzer Space Telescope (Werner et al. 2004). The images are 
publicly available and are downloaded from the Spitzer Space 
Telescope Archive using the Leopard package 2 . These ob- 
servations (dataset archive ID: AOR 12485 120) were obtained 
within the GTO Science Program 125 (PI: G. Fazio) on 28 
November 2004. We retrieved the Post Basic Calibrated Data 
(pbcd) for all four IRAC channels (3.6 /im, 4.5 /an, 5.8 
fjm and 8.0 /im), which are already fully co-added and cal- 
ibrated. After visual inspection of the images we carried 
out PSF-photometry with the daophot package provided 
within the I RAF 3 environment. In the case of the region of 
NGC 602/N 90, PSF photometry is required as the field is too 
crowded for aperture photometry to be appropriately applied. 
Consequently, for each filter we constructed a reference PSF 

2 Leopard can be downloaded from 

http://ssc. spitzer. caltech. edu/propkit/spot/ 

3 I RAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation, (http : / /iraf . noao . edu / 1 



by combining the PSF of at least 5 objects with a high signal 
to noise ratio and we used this PSF to fit all objects identified 
with a 3-er confidence level over the local background. 

We converted the pixel values in each filter from MJy sr -1 
to DN s" 1 according to "IRAC Data Handbook" version 5.1.1 
and we computed the corresponding magnitudes as m = 
-2.5 log(/) + Azp with / denoting the flux measured in 
DNs" 1 and Azp being the zero point for each filter. Zero 
points were taken from Hartmann et al. (2005): 19.66 (3.6 
/im), 18.94 (4.5 /mi), 16.88 (5.8 /im) and 17.39 (8.0 /im). 
For the initial PSF fit we used a PSF size of 2 pixel, and 
applied aperture corrections as described in the "IRAC Data 
Handbook" to obtain the final magnitudes. Finally, we cross- 
matched the obtained catalogs for each filter to identify the 
sources that are detected in more than one IRAC bands, by 
allowing a maximum offset of 2 pixels for the center of each 
point source. 

3. DESCRIPTION OF THE MID-INFRARED OBSERVATIONS OF 
NGC 602/N 90 

Spitzer covers a suitable wavelength range (3.6 - 160 /im) 
for defining the dust properties of the interstellar medium 
(ISM). The mid-infrared channels of IRAC (3.6, 4.5, 5.8 and 
8 /im) are good tracers of the smallest grains through Poly- 
cyclic Aromatic Hydrocarbon (PAH) emission. Very small 
grains (typical sizes 12 - 150 A) and PAHs (typical sizes 4 - 
12 A; Desert et al. 1990) are very efficient in heating the gas 
through the photoelectric effect (Bakes & Tielens 1994) and 
therefore their analysis is very important for understanding 
the thermodynamics of the ISM. 

IRAS data on both LMC and SMC suggest a deficit of the 
smallest dust grains in the low-metallicity environments of 
these galaxies due to the lower 12 /zm emission in compar- 
ison to the Milky Way and other higher metallicity galaxies 
(Schwering 1989; Sauvage et al. 1990). A study in the SMC 
by Stanimirovic et al. (2000) concluded that the dust in the 
SMC is dominated by large grains. The lower abundance of 
PAHs and very small grains in the SMC is also supported by 
detailed modeling of SMC UV absorption and FIR emission 
data by Weingartner & Draine (2001), Li & Draine (2002), 
and Clayton et al. (2003). A different grain size distribution 
with respect to the Galaxy, or grain destruction, possibly due 
to the intense UV radiation and/or by supernova shocks have 
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FIG. 2. — Left: IRAC CMD of template/model photometry of Cohen (1993) and Whitney et al. (2004) for IR sources in the Milky Way. Symbols code different 
types of sources throughout the SMC: Main-sequence O stars, RSGs, O-rich and C-rich AGB stars, H II regions, and YSOs. Right: The corresponding CMD of 
sources found in the region of NGC 602/N 90 with our IRAC photometry. The comparison of these CMDs denotes that there are 27 objects in our photometry 
(located at the right part of the CMD on the right with colors [3.6] — [8.0] > 2 mag) that match the type of galactic YSOs. 



been proposed to explain the very small grain paucity in both 
galaxies. The absence of PAHs and very small grains has a 
profound influence on the gas heating and the existence of 
cold and warm phases in the interstellar medium (Wolfire et 
al. 1995). 

The images of NGC 602/N 90 taken with Spitzer, which we 
present here, can be used for a qualitative analysis on the dust 
properties and its relation to stellar sources of radiation in the 
region of NGC 602. A color composite image from IRAC 
channels 1, 2 and 4 is shown in Fig. \T\ The 3.6 and 4.5 /im 
wave-bands are mostly sensitive to stellar photospheres and 
very hot circumstellar dust, and therefore the blue (3.6 /im) 
and green colors (4.5 /im) in Fig. Q]reveal the stars, whereas 
the red color (8 /im) shows mainly the diffuse dust emission. 
Specifically, the origin of diffuse emission at 3.6 /im is likely 
to originate partly from bound-free transitions in the ionized 
gas with some contribution from the 3.29 /im emission near 
molecular clouds (thought to be due to C-H bond stretching 
in PAHs) and very small dust grain emission (e.g. Engelbracht 
et al. 2006). The 4.5 /im diffuse emission is a combination 
of bound-free continuum and Brackett-a recombination emis- 
sion in H II regions. Both the 5.8 and 8.0 /im bands show 
molecular material in the interstellar medium (ISM) and cir- 
cumstellar envelopes, as well as increasingly faint stellar pho- 
tospheres. The origin of the 5.8 /im diffuse emission is most 
likely very warm dust emission (T ~ 600 K). It encompass the 
5.6 and 6.2 /im PAH band emission. The origin of the 8 /*m 
diffuse emission is most likely the very bright emission com- 
plex at 7 - 9 /im thought to be dominated by C-C stretching 
modes (with some contribution from in-plane C-H bending) 
of the bonds in PAHs and very small carbonaceous dust grains 
(Jager et al. 2006). 

4. POINT-SOURCES DETECTED WITH SPITZER IN NGC 602/N 90 

Our photometry revealed 57 objects detected in all four 
IRAC channels. 44 sources were identified in channels 1, 2 
and 3, but not in channel 4, while only two in channels 2,3 and 
4 and not in channel 1 . For the classification of the sources we 
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FIG. 3. — Spitzer/IRAC Color-color diagram of the sources detected in 
our photometry in the region of NGC 602/N 90. The vertical dashed line 
corresponds to the limit [3.6]-[8.0]=1.2 mag, set by Bolatto et al. (2007) 
as a criterion for the selection of candidate YSOs similar to those modeled 
by Whitney et al. (2003). We refine our selection of candidate YSOs in the 
region to 22 point-sources, which are redder than this limit. 



make use of the [3.6]-[8.0], [8.0] color-magnitude diagram 
(CMD) in comparison with templates and models compiled 
for a wide range of IR sources in the Milky Way. These tem- 
plates have been developed by Cohen (1993) for IRAS data 
on galactic IR sources, adapted for the IRAC and 2MASS 
color classification scheme for the Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire (GLIMPSE) Spitzer project 
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FIG. 4. — Color composite image of separate exposures made by ACS/WFC on Hubble Space Telescope in two broad-band filters, F555W (V) and F814W (/), 
and one narrow-band filter, F65SN (Ha + [N II]) of the region NGC 602/N 90. North is up and east to the left. The data of these observations are the subject 
of an investigation of the young stellar populations in the region (Schmalzl et al. 2007). Red boxes indicate the positions of the candidate YSOs detected with 
our IRAC photometry. This picture clearly indicates that current star formation takes place at the periphery of the central association NGC 602. Probably, star 
formation started at the association (where the highest concentration of PMS stars has been found; Schmalzl et al. 2007) and propagated outwards, mainly to the 
east and south, with events of ongoing clustered star formation still taking place along the dust ridges. 



(Benjamin et al. 2003). We also make use of the YSO models 
that have been developed by Whitney et al. (2004) also for 
GLIMPSE. Both papers include several types of IR sources 
such as main-sequence stars, red giants, O-rich and C-rich 
AGB stars, H II regions, and YSOs. In Fig. |2](/e/f), we show 
the position of such sources in the [3.6]-[8.0], [8.0] CMD 
according to the templates and models. We also show the cor- 
responding CMD for the point-sources we found in all four 
IRAC channels in the region of NGC 602/N 90 (Fig. [2 right). 
From the comparison of the CMDs of Fig. |2]it is shown that 
there are 27 objects in our photometry with positions in the 
[3.6]-[8.0], [8.0] CMD similar to galactic YSOs of class 0, 1 
and II. However, Class YSOs are very rare due to their short 
evolutionary phase, and such objects would be extremely faint 
in the IRAC bands at the distance of the SMC. Consequently, 
one should consider that, although few of our sources match 
the positions of Class YSOs in the CMD according to the 
models, probably there are no YSOs of this class in our sam- 
ple. 

A criterion adopted by Bolatto et al. (2007) within the 
S 3 MC Survey (see also §1) for the selection of YSOs in their 
sample is based on very red [5.8]-[8.0] color index (> 1.2 
mag) of bright 5.8 /im sources. Support for their selection 
is lent by the location of these sources in the [3.6] - [8.0], 
[8.0] CMD, and the color-color diagrams in near-, mid- and 
far-infrared wavelengths as modeled by Whitney et al. (2003) 
for different objects. Bolatto et al. note that the highly red- 
dened 8.0 /im sources occupy the region in those plots that 
is predicted for late Class to Class II YSOs with different 
geometrical parameters. We further select the most proba- 
ble YSOs in the region of NGC 602/N 90 by applying this 
criterion to the sources detected with our photometry. The 
[5.8]-[8.0], [3.6]-[4.5] color-color diagram of our sources is 
shown in Fig. [3]with the limit of Bolatto et al. (2007) for can- 
didate YSOs overplotted as a vertical dashed line. There are 
24 sources in our photometry that have colors [5.8]-[8.0] > 
1.2 mag, meeting this selection criterion. All, except two of 
them, have been selected as candidate YSOs also according to 
their positions in the [3.6]-[8.0], [8.0] CMD (Fig.©. Conse- 
quently, we refine the number of candidate YSOs in the region 



of NGC 602/N 90 to 22 sources, which meet both criteria. 

The classification scheme of Bolatto et al. predicts only ~ 
2% contamination of this part of the mid-infrared color-color 
diagram by B stars with 24 /im emission (their Figure 12). 
Other types of evolved sources, such as carbon stars, RSGs 
and O-rich AGB stars, are located to bluer colors. In addition, 
Bolatto et al. estimated the contamination by background, un- 
resolved galaxies, which can be seen throughout the SMC at 
infrared wavelengths. They expect 20 contaminants in their 
sample of 282 candidate YSOs (~ 7%), which translates to 
1.5 sources in our sample. However, detailed inspection of 
the candidate YSOs on the ACS image of NGC 602/N 90 re- 
vealed no suspicious background galaxies in our sample. The 
characteristics of the candidate YSOs revealed with our pho- 
tometry in the region of NGC 602/N 90 are presented in Table 
[U where we give the ID number (from our IRAC photometric 
catalog), the position and the magnitudes and corresponding 
photometric uncertainties in the four IRAC channels for each 
candidate YSO. 

It should be noted that the classification of candidate YSOs 
cannot, in any case, be reduced to a single criterion, and there- 
fore the technique described above should be considered a 
first-order selection. Moreover, taking into account that the 
selection of objects with [5.8]-[8.0] > 1.2 mag may compro- 
mise the actual number of YSOs, this limit should be consid- 
ered a tentative one. Indeed, Allen et al. (2004) plot known 
stars with YSO classifications in the IRAC color-color dia- 
gram (their Fig. 4) and their Class I YSOs clearly stretch 
to red [3.6]-[4.5] colors with [5.8]-[8.0] colors less than 1.2 
mag. Our selection, based on Bolatto et al. (2007) scheme, 
apparently causes such sources to be left out, but its use is 
being considered in order to reduce as much as possible the 
number of misidentifications. Furthermore, sources with col- 
ors [3.6]-[4.5] w 0.0 mag and [5.8]-[8.0] > 1.2 mag are not 
well represented in Whitney et al. (2003, 2004) nor in Allen et 
al. (2004). Based on their colors, some of them are suspected 
to be small H II regions, but there are no high resolution radio 
observations of N 90 (unlike N 159 in the LMC; Jones et al. 
2005), which would provide information on the ionizing flux 
for these sources and help us clarify their nature. 
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FIG. 5. — Color composite images constructed from our ACS observations of regions 18" X 18" wide centered on the loci where candidate YSOs have been 
detected. Characteristic examples of 13 candidate YSOs from TablefT]have been selected. The drawn circles, corresponding to radius 1"2 indicate the positions 
of the candidate YSOs. Blue circles denote candidate YSOs for which optical sources have been found to coincide with, while red circles the ones for which no 
optical sources have been identified. 



TABLE 1 

Candidate YSOs detected in all four IRAC channels in the region of NGC 602/N 90. The numbering system of the YSOs is from 

OUR PHOTOMETRIC IRAC CATALOG. CELESTIAL COORDINATES ARE GIVEN IN J2000. THE MAGNITUDES OF EACH CANDIDATE YSO IN ALL CHANNELS 
WITH THE CORRESPONDING ERRORS ARE GIVEN IN COLUMNS 4 THROUGH 1 1 . THE LAST FIVE COLUMNS REFER TO THE OPTICAL SOURCES FOUND WITH 
ACS TO FORM POSSIBLE COMPACT CLUSTERS WITHIN A RADIAL DISTANCE OF 1 "2 AROUND EACH OF THE CANDIDATE YSOS. THE NUMBERS OF THE 
DETECTED OPTICAL SOURCES, THE MAGNITUDES OF THE BRIGHTEST SOURCE (F555W bri , F814W bli ), AS WELL AS THE CUMULATIVE MAGNITUDES OF 
ALL STARS AROUND EACH CANDIDATE YSO (F555W to t, F814Wiot) ARE GIVEN IN COLUMNS 13, 14, 15 AND 16 RESPECTIVELY. 



ID 




RA 


DEC 


[3.6] 


± 


[4.5] 


± 


[5.8] 


± 


[8.0] 


± 


No 


F555Wbri 


F814W bri 


F555W tot 


F8I4W, ot 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


1 


01 


29:35.66 


-73:33:39.61 


15.90 


0.22 


16.24 


0.30 


13.64 


0.24 


11.99 


0.15 


2 


24.62 


22.85 


24.46 


22.63 


2 


01 


29:25.95 


-73:34:46.69 


16.52 


0.23 


16.82 


0.32 


14.17 


0.31 


12.50 


0.19 


5 


25.45 


23.62 


24.66 


22.45 


5 


01 


29:26.52 


-73:32:30.31 


16.27 


0.15 


16.45 


0.21 


14.57 


0.29 


12.91 


0.15 


2 


17.84 


17.72 


17.82 


17.66 


22 


01 


29:44.77 


-73:33:25.31 


16.46 


0.21 


16.50 


0.26 


14.70 


0.32 


13.06 


0.17 


15 


23.26 


21.76 


21.82 


20.28 


24 


01 


29:39.23 


-73:33:18.80 


16.66 


0.20 


16.69 


0.31 


14.45 


0.36 


12.79 


0.22 


5 


24.54 


21.85 


23.12 


20.93 


26 


01 


29:19.80 


-73:33:12.46 


16.06 


0.23 


16.04 


0.29 


13.84 


0.30 


12.17 


0.17 












29 


01 


29:35.55 


-73:33:30.20 


14.13 


0.09 


14.09 


0.11 


12.22 


0.13 


10.61 


0.08 


8 


17.49 


17.30 


16.73 


16.65 


34 


01 


29:23.08 


-73:32:51.51 


16.00 


0.16 


15.93 


0.17 


14.28 


0.26 


12.88 


0.13 


4 


22.82 


20.95 


22.06 


20.13 


35 


01 


29:22.34 


-73:32:48.17 


16.67 


0.27 


16.60 


0.36 


14.63 


0.35 


12.95 


0.22 












36 


01 


29:28.61 


-73:34:12.14 


16.88 


0.31 


16.80 


0.42 


15.32 


0.50 


13.68 


0.29 


1 


25.15 


24.08 






38 


01 


29:34.08 


-73:34:11.66 


16.67 


0.28 


16.58 


0.39 


14.84 


0.39 


13.10 


0.18 












42 


01 


28:47.20 


-73:31:31.69 


16.70 


0.24 


16.54 


0.32 


15.37 


0.51 


14.02 


0.29 












43 


01 


29:41.38 


-73:36:24.53 


15.92 


0.11 


15.71 


0.15 


14.67 


0.33 


12.87 


0.11 












45 


01 


29:43.34 


-73:33:43.21 


16.50 


0.20 


16.21 


0.20 


14.19 


0.29 


12.97 


0.18 


2 


21.90 


21.22 


21.89 


21.16 


46 


01 


29:30.92 


-73:34:05.74 


16.51 


0.23 


16.17 


0.27 


14.80 


0.38 


13.36 


0.21 


1 


22.45 


20.92 






49 


01 


29:20.44 


-73:33:16.90 


16.68 


0.21 


16.19 


0.23 


15.16 


0.50 


13.65 


0.21 


1 


23.16 


21.55 






50 


01 


29:37.52 


-73:33:51.34 


15.00 


0.18 


14.50 


0.20 


13.09 


0.22 


11.70 


0.19 


1 


20.96 


20.42 






51 


01 


29:06.45 


-73:33:48.73 


13.91 


0.08 


13.35 


0.10 


12.15 


0.11 


10.69 


0.06 


10 


20.16 


19.56 


20.05 


19.19 


52 


01 


29:30.21 
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NOTE. — Objects No 54 and 57 being very close to each other, have been identified as distinct sources only in the 3.6 /^m band. Our photometry could not resolve them in IRAC 
channels 2,3 and 4, and therefore they appear with the same corresponding magnitudes. Inspection of their ACS image (Fig. [5} shows indeed the existence of multiple objects very close 
to each other. Around each of the candidate YSOs 54 and 57 different PMS stars have been found with only four common stars. Objects No 42 and 43 are located outside the field 
observed by ACS, and consequently it was not possible to identify any optical sources for them. 



4.1. Optical Counterparts of the Candidate YSOs 

We use the results of the photometry by Schmalzl et al. 
(2007) on ACS/WFC imaging of NGC 602 to identify the op- 
tical counterparts of the candidate YSOs found with Spitzer 
located in this region. We cross-correlate the positions of the 
candidate YSOs with the ones of all stars found with ACS. 
For the comparison of the two catalogs we take into account 
the fact that the resolution provided by IRAC is almost 25 



times lower than the one of WFC, and therefore the selected 
area around each candidate YSO to be searched for optical 
counterparts in most of the cases revealed more than one op- 
tical sources. Specifically, we selected a search circle with 
diameter of 2. "4, which corresponds to the size of two IRAC 
pixels, a typical FWHM of the point-sources found with our 
PSF photometry. In TableQ]the results of this search are given. 
Specifically, in the last five columns of this table we give the 
number of stars from the ACS photometry, found within a ra- 



D. A. Gouliermis, S. P. Quanz, & Th. Henning 



' 6 




1 2 

V— I [mag] 



FIG. 6. — V — I, V Color-Magnitude Diagram of all the optical sources 
from our ACS photometry, the loci of which are found to coincide (within a 
radial distance 1"2) with the location of the candidate YSOs identified with 
our IRAC photometry in the region of NGC 602/N 90. Evolutionary tracks 
for PMS stars from Palla & Stahler (1999) for masses 0.2, 0.4, 0.6, 0.8, 1.0, 
1.2, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 5.0 and 6.0 M are overplotted. This diagram 
shows that these optical sources are PMS stars, which cover a wide range of 
masses. PMS isochrone models for 0.1 and 5 Myr and the ZAMS are also 
plotted with thick solid lines. These models limit the age of the PMS stars, 
and consequently the timescale for the corresponding star formation events, 
to ^ 5 Myr. 



dial distance of 1 ."2 (2 IRAC pixels) around the center of each 
candidate YSO. We also give the magnitude of the brightest 
optical source in the filters F555W (V) and F8 14W (I) respec- 
tively, as well as the corresponding total magnitudes from all 
optical sources found within the searched area around each 
candidate YSO. In TableQ]there are five candidate YSOs for 
which no optical counterpart was found. However, only three 
of them (No 26, 35 and 38) are located within the ACS/WFC 
field-of-view for which the cross-corelation was possible. The 
remaining two candidate YSOs (No 42 and 43) are located 
outside this field-of-view, since the selected area observed by 
IRAC (Fig. [TJ is a bit larger than the one covered by ACS. 

5. CURRENT STAR FORMATION IN THE AREA OF NGC 602/N 90 

5.1. Clustered Star Formation Events 

As shown in Table [T] almost all of the detected candidate 
YSOs in the area of NGC 602/N 90 are related to stars visible 
at optical wavelengths. In half of the cases there are more than 
two optical sources related to each candidate YSO (within a 
region corresponding to two IRAC pixels). The loci of the 
detected candidate YSOs are shown marked with red squares 
on a color composite image from the observations taken with 
ACS/WFC of the area of NGC 602/N 90 in Fig. g] In this fig- 
ure it is shown that there is a higher concentration of candidate 
YSOs to the east and south, where also the nebula is brighter 
and where the brighter 8/im emission is found in the Spitzer 
images (Fig. [TJ. More detailed color composite images from 



ACS imaging of the loci around 13 selected candidate YSOs 
are shown in Fig. [5] Circles of radius corresponding to 1"2 
mark the centers of the candidate YSOs. In this figure it can 
be seen that there are candidate YSOs coinciding with stellar 
clusterings, which could be true compact clusters (objects No 
22, 24, 29, 54/57). The nature of these compact clusters and 
if they are signatures of on-going clustered star formation are 
two very important issues. 

We attempt to address these issues by plotting the CMD of 
all stars from our ACS photometry, that are associated with 
the candidate YSOs found with IRAC. The part of a V-I, V 
CMD, which is expected to be occupied by pre-main sequence 
(PMS) stars of star-forming regions in the MCs observed with 
Hubble is fairly known (e.g. Gouliermis et al. 2006, Nota et 
al. 2006). This CMD for the optical sources related to the can- 
didate YSOs in the region of NGC 602/N 90, plotted in Fig. 
|6l shows that all of the stars found to coincide with the candi- 
date YSOs in our sample are actually PMS stars, suggesting 
that either we probably observe the existence of compact PMS 
star clusters recently formed (or currently under formation) 
around typical YSOs, or that the sources identified as candi- 
date YSOs are themselves unresolved young PMS clusters. In 
either case it is almost certain that these stellar concentrations 
demonstrate possible events of extragalactic clustered star for- 
mation on the act. If so, it would be interesting to know how 
long ago this formation took place. 

To answer this question, on the CMD of Fig. [6] we also plot 
PMS evolutionary tracks from the models by Palla & Stahler 
(1999), and two indicative isochrones (for 0. 1 and 5 Myr) and 
the Zero Age Main Sequence (ZAMS) from the same models. 
Simulations of the positions at the V-I, V CMD of all PMS 
stars found in the region of NGC 602/N 90 by Schmalzl et 
al. (2007) indicate that reddening, variability and binarity are 
important factors, which cause a spread of the loci of PMS 
stars in this CMD and that even modest interstellar redden- 
ing (Schmalzl et al. find E(B-V) ~ 0.04 - 0.06 mag for the 
association NGC 602) can make these stars appear younger 
(being "pushed" to redder colors). Considering these results 
and the plot of Fig. |6l one may conclude that the PMS stars, 
which coincide with the candidate YSOs of our sample do not 
seem to be older than about 5 Myr, giving an upper limit to 
the time-scale within which these clusters were formed. 

For a more detailed analysis one should concentrate on the 
stellar content of each of the individual clusters, but in most 
of the cases very low numbers of stars per cluster do not allow 
such an analysis. We selected the most populous among the 
identified clusters (the ones corresponding to candidate YSOs 
No 22, 29, 51 and 54) to plot their individual CMDs, but the 
stellar numbers were still not enough to define a more narrow 
age range. PMS stars are mainly red sources, and therefore 
near-infrared imaging with cameras equipped with Adaptive 
Optics (AO) systems would provide the necessary resolution 
for a richer sample of PMS stellar members in these clusters. 
Indeed, high resolution near-infrared observations of these 
clusters would certainly be very useful, because they would 

1) provide better statistics with higher numbers of detected 
members, enough to plot the CMDs of individual clusters and 

2) fill the gap between optical and mid-infrared wavelengths, 
providing near-infrared magnitudes for these PMS stars. The 
later information will offer the opportunity to construct accu- 
rate spectral energy distributions (SEDs) of candidate YSOs, 
to be compared to models. 

As far as candidate YSOs that coincide with no more than 
two PMS stars concerns, in Fig. [5] we also show several dif- 
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ferent cases of "single" candidate YSOs (e.g. objects No 1, 5, 
49, 53), but as their images show it is unclear if such sources 
are indeed single isolated objects or embedded clusters (as the 
faint magnitudes given in Table Q] suggest), or even one mas- 
sive embedded protostar. Modeling of the SEDs of such ob- 
jects is quite difficult, in the sense that many different models 
seem to fit the SEDs, due to their small number of points, and 
near-infrared AO observations would certainly add three more 
points to the SEDs, throwing more light toward this direction 
as well. It is worth noting that Fig. [5] also demonstrates few 
interesting cases of candidate YSOs (objects No 26, 49, 34, 
35, 46 and 53), both with and without optical counterparts, 
which are located at the edge of features similar to the "Pil- 
lars of Creation" or "Elephant Trunks" in the galactic H II 
region M 16 (the Eagle Nebula; Hester et al. 1996). 

5.2. Star-Formation Scenario for the Region of 
NGC 602/N 90 

The sample of candidate YSOs derived from our IRAC pho- 
tometry in the region of NGC 602/N 90 (TableHJ is based on a 
single classification scheme. As a consequence, although it is 
almost certain that this sample does include true YSOs, there 
is always a possibility that it also includes sources misiden- 
tified as YSOs, like small H II regions, or that it is quite in- 
complete not including all possible true YSOs, like the Class 
I YSOs in the sample of Allen et al. (2004). The identifi- 
cation of PMS stars in the areas of the candidate YSOs as 
their optical counterparts indicates that half of these objects 
are not single sources. But, considering the difficulty of mod- 
eling the SEDs of these objects, it is not straightforward to 
clarify which of them represent embedded proto-clusters. In 
the case of candidate YSOs coinciding with one or two PMS 
stars, there is no information if there should be more opti- 
cal counterparts, which simply were not detected in our ACS 
photometry, or if these are indeed single proto-stellar objects. 

In any case, the composite ACS image of the area of 
NGC 602/N 90 shown in Fig. H] provides evidence that se- 
quential star formation does take place in the area and espe- 
cially on the dust ridges of a cavity, which seems to extend 
outwards. This feature is most probably the product of pho- 
toionization by the young association, NGC 602, located al- 
most at its center. The combination of both ACS (Fig. HI and 
IRAC (Fig. [TJ observations provides a more complete image 
of the star formation in the area. Specifically, taking into ac- 
count the loci of the sources classified as candidate YSOs, the 
ISM features shown in the image from ACS, and the 8 /im 
emission shown in the IRAC image, one may conclude that in 
the observed area the stellar association NGC 602 ionizes the 
surrounding material giving birth to the H II region, N 90, and 
producing a cavity in the remaining of its parental molecular 
cloud. This event triggers sequential star formation not later 
than ~ 5 Myr ago, which seems to propagate mostly to the 
east and south. This second star formation event may be char- 
acterized by higher a concentration of YSOs, most of which 
may be related to compact PMS clusters. This scenario is also 
supported by the study of the PMS population in NGC 602 
and its surroundings by Schmalzl et al. (2007). 

6. CONCLUSIONS 

We present our results from Spitzer/IRAC photometry on 
the star-forming region NGC 602/N 90 in the Small Magel- 
lanic Cloud. 

We performed PSF photometry on IRAC images taken in 
four mid-infrared channels (3.6 /im, 4.5 /im, 5.8 /im and 8.0 



/im) centered on the H II region N 90. We identified 57 ob- 
jects in all four channels and 44 in channels 1, 2 and 3. We 
classified the sources found in all four IRAC channels from 
their loci in the [3.6]-[8.0], [8.0] color-magnitude diagram, 
based on templates and models for different red sources in 
the Milky Way, and the [5.8]-[8.0], [3.6]-[4.5] color-color 
diagram, according to the classification scheme developed by 
Bolatto et al. (2007) for SMC YSOs. 

According to our classification there are 22 candidate YSOs 
in the region of NGC 602/N 90. We use our results from PSF 
photometry in the filters V and / of the same region from 
imaging with the Wide-Field Channel of ACS on-board the 
Hubble Space Telescope (Schmalzl et al. 2007) to identify 
optical sources, which coincide with the candidate YSOs. All 
these optical sources are identified as PMS stars. We found 
that almost half of the candidate YSOs match with no more 
than two PMS stars, but examination of their loci on the ACS 
images does not exclude the possibility that there may be 
more highly embedded sources or not well resolved also in the 
optical. The remaining candidate YSOs coincide with more 
than 3 PMS stars, and four of them seem to include small 
compact PMS clusters of 7 (or more) stars. We found that 
these stars are probably formed earlier than 5 Myr ago, ac- 
cording to PMS evolutionary models. We could not find any 
optical counterpart for five YSOs, two of which are located 
outside the field observed with ACS. 

From the positions of the detected candidate YSOs, we con- 
clude that they are probably the products of a sequential star 
formation process, which presumably propagates mostly to 
the east and south, on the ridges of the molecular cloud, as 
they formed from the photo-ionizing process of the stellar as- 
sociation NGC 602 located in the center of the apparent cavity 
of the cloud. 

Finally, it should be mentioned that our sample of candi- 
date YSOs may be contaminated by other objects like small 
H II regions, or it may not include some true YSOs. Based on 
our results we stress the importance of high-resolution near- 
infrared imaging of the candidate YSOs in order 1) to resolve 
more faint red stellar sources and detect a higher number of 
related PMS stars, and 2) to fill the wavelength gap between 
optical and mid-infrared by providing near-infrared measure- 
ments of all detected sources, so that their SEDs will be thor- 
oughly modeled, as well as of high-resolution radio observa- 
tions, which will provide additional information on the ioniz- 
ing flux for the identified objects. 
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